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pH ■ Abstract. At low Reynolds numbers, the hydrodynamic interaction between 

dumbbells driven by an external rotating field can be attractive or repulsive. Dumbbells 
^ ' of dissimilar asymmetric shape or different coupling to the external field undergo 

^ ■ conformational rearrangements that break the time reversal symmetry. The parameter 

(^_| I ranges leading to attraction or repulsion are explored numerically. The results of 

' our simulations suggest that rotating fields may be a useful avenue for the assembly, 

disassembly, and sorting of particles of different shape as well as for the study of 
^ ' collective micro-swimmers. 

1. Introduction 

> 

^ I Hydrodynamic interactions in suspensions of active or externally driven nanoparticles 
^ ' or microorEranisms are involved in collective, biological, chemical and technical motion 
^1 [Il[2l[3llll[5l[6l[7l[8l[9]of fascinating complexity. The propulsion of micro- 

O I particles and microorganisms differs substantially from the motion of macroscopic 
^ ' objects like fishes or birds due to the small scale motion at low Reynolds numbers 
[2]. As a consequence swimming macro- and microorganisms transform chemical into 
^ ■ translational or rotational energy in different ways. Considerable progress has been 
' made in understanding the propulsion of microorganisms by studying artificial active 

swimmers [91110,111] or active rotors [I2]. Self-organization phenomena of assemblies of 
micro-swimmers can teach us how hydrodynamic interactions synchronize the motion of 
swimmers or cilia The importance of hydrodynamic interactions (HI) at small Reynolds 
numbers has been emphasized by showing that they may cause a lift force during the 
motion of vesicles close to a wall [131 El IIS]- They also lead to cross-streamline 
migration of droplets, vesicles and deformable bead-spring models in Poiseuille flow 
[T6 l [T7 1 [T8]. Additionally, the interplay between Brownian motion of small particles and 
hydrodynamic interactions can cause surprising anti-correlations and cross-correlations 
[T9l[20]. Brownian dynamics of hydrodynamically interacting polymers may cause elastic 
turbulence that is used for mixing in microfiuidic devices [211 122] • 
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Macroscopic rotors experience the attractive Magnus effect that can easily be 
understood using the concept of dynamic pressure. The attraction or repulsion between 
rotors in the limit of Stokes flow is less obvious. Magnetically driven rotors can attract 
each other due to their magnetic interaction [6l [23l [2l] . They form dynamically induced 
self-organized structures depending on the particles' contours Magnetic dipole- 
hexapol interactions may be used for separating rotors with respect to their shape 
[25] . Rotors also play an important role for bacterial motion such as the propulsion 
of Escherichia coli [3] with its rotating flagella. The dynamics of single fiagella has 
recently been modeled in numerical calculations [26] . There hydrodynamic interactions 
are the source of the propulsion the fiagella create. One fundamental problem for active 
and driven low Reynolds number motion is the understanding of the consequences of 
hydrodynamic interactions on the collective motion of an assembly of active components. 

The aim of this work is to explore the effect of hydrodynamic interactions on 
two rotating dumbbells. Symmetry considerations forbid a net conformational change 
over one rotational period for similarly shaped dumbbells with similar couplings to the 
external field. One major question explored here is the influence of differing shapes 
on the hydrodynamic interactions between two particles. Therefore we investigate two 
interacting dumbbells of dissimilar shape and coupling to an external rotating field. 

We demonstrate by numerical calculations that for dissimilar rotating dumbbells 
the time reversal symmetry is broken. As a consequence both dumbbells can attract or 
repel each other. Transitions between dumbbell attraction and repulsion are induced 
by changing the torques or the dumbbell shapes. Phase diagrams separating regimes of 
attraction from regions of repulsion between the dumbbells are constructed for generic 
sets of parameters. 

We suggest to test our theoretical exploration with experiments on anisotropic 
birefringent particles. Such particles may be rotated by circularly polarized light [271 [28]. 
Another possibility to apply a torque is to place paramagnetic or ferromagnetic particles 
in a rotating magnetic field. Like in nuclear magnetic resonance, magic angle spinning 
may be used to eliminate time averaged dipole-dipole interactions that might overshadow 
the hydrodynamic attraction or repulsion. 

2. Model 

We investigate the motion of a single dumbbell and of two hydrodynamically interacting 
dumbbells in a fluid of viscosity rj in the low Reynolds number limit. Due to a torque 
created by a rotating external field / the dummbbells perform a planar rotational 
motion. Each dumbbell is composed of two beads which are connected by a bar of 
constant length as indicated in figure [H 

The beads with the effective hydrodynamic radii a-i {i = 1,2,3,4) at the positions 
Ti {i = 1,2,3,4) are described as point particles having drag coefficients Q = Gnrjai. 
The vectors r'21 = r2 — vi and r^^ = — describe the two dumbbell axes, which 
we assume to be of constant length b = \r2i\ = |?"43|. The hydrodynamic centers of the 
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dumbbells, C21 and C43, are given by 

(32 



C2I = Ti + 



ai + 02 



r2i 



and 



C43 = ^"3 + 



04 



as + 04 



(1) 





Figure 1. Sketch of two dumbbells with unequal asymmetries, ai ^ and 02 = 04. 
The dumbbells are rotated by a driving field f{t), cf. Eq. ([5]). is the angle between 
the axes r2i and of the dumbbells. The crosses indicate the hydrodynamic centers 
C21 and C43 of the two dumbbells. The distance between them is denoted by Ch- 



The equation of the over-damped motion is given by 
r-i = HijFj , 

with the Rotne-Prager mobility matrices for unequal spheres [29l [30] 

^ -I for i = j, 

I + ( 1 - ) rijr,j for z 7^ j. 



(2) 



H, 



1 



1 



3^ 



(3) 



Tij/rij and the external forces -Fj (j 



the connection unit vectors f-ij = rij/\rij\ 
1, 2, 3, 4) acting on the beads. In our model the rotation of the dumbbells is caused by 
the potential 



V = -V2lf2lf - Kt3^43/, 

with the vector 

'cos(a;t)~ 
sin(a;t) 

rotating at a frequency u in the xy plane. The forces Fi 
expressed in terms of the rotating field / as 



(4) 



(5) 



-ViV in Eq. (E]) can be 



Fi = — [(fsi X /) X fai] and F3 = — [(^^43 x /) x 7^43] . (6) 

^"21 ^"43 

The potential in Eq. yields the following asymmetric relations between the forces 
Fi = -ViV = V2V = -F2 and F3 = -V3I/ = V4V = -F4. According to this 
property together with Eq. the dumbbells are force free. Additionally the forces 
acting on the individual beads are perpendicular to the dumbbell axes, which results in 
a rotation of the dumbbells. 

The form of the potential V for the driving field is similar to that of magnetic dipoles 
in an external magnetic field where the interaction between the dipoles is neglected. 



v. 



43 
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But other means to propel the dumbbells are possible. For example, non-spherical 
birefringent particles can be rotated by circularly polarized light. The magnitude of the 
torques imposed on the dumbbells can easily be tuned by the applied power of the light 
[271 [28]. 



3. Dumbbell dynamics 

The balance between the driving and the viscous torques leads to phase angles (j)2i 
and 043 between the dumbbell orientations Vij and the driving field f(t) with cos(0jj) = 
Tij- f /{vij- f\. The dynamics of the dumbbells is analyzed in this section. The conditions 
under which both dumbbells attract or repel are discussed. 



3.1. Single dumbbell 

The two beads of a single dumbbell (V43 = and V21 = Vq) move on circular trajectories 
around their hydrodynamic center C21. The radii Ri^2 of those trajectories depend on 
the effective hydrodynamic radii ai^2 of the two beads and the length b of the dumbbell 
axis. In the Oseen approximation they are given by 

R, = and i?2 = . (7) 

ai + a2 ai + a2 

For a sufficiently strong external field the dumbbell axis r2i rotates synchronously 
with the driving field f{t) at the same frequency Q = u after a transient period. The 
stationary phase angle 02i is given by the expression 

Tcrjub'^ 



sm (p2i 



Vo 



1 ^ 1 1 / -j^ , + ^2 



(8) 



6ai Qa2 4h \ 36^ 

021 tends to zero for decreasing values of the ratio w/Vq, i. e. either for a decreasing 
frequency or for increasing values of the coupling Vq. If the expression on the right 
hand side in Eq. ([8|) takes values outside of the interval [—1, 1], the driving torque is too 
small to enforce a synchronous rotation. These analytical results are in agreement with 
numerical integrations of Eq. ([2]). 



3.2. Two dumbbells 

In the case of two rotated dumbbells each one experiences the perturbed liquid fiow 
created by the respective other one. Due to this interaction the hydrodynamic centers 
Cij are set in motion. In the mean this results in a circular or a spiral like motion as 
indicated in figure [2J In the following considerations we fix the effective radii of one 
bead of each dumbbell at the value 02 = 04 = 1 and the length of the axes of both 
dumbbells at 6 = 3 for simplicity. The (asymmetric) shapes of the dumbbells and the 
couplings to the external field can be adjusted independently by varying the bead radii 
ai and 03 as well as the coupling parameters V21 and V43. The viscosity rj, the angular 
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frequency u, the coupling Vij and the length b form a dimensionless group rjub^/Vij. It 
is therefore sufficient to vary the coupling and keep the viscosity rj = 10 and the angular 
frequency u = 5 x 10~^ fixed. 

The complete conformation of the system (up to a rotation of the whole system) can 
be described by the distance D{t) between the hydrodynamic centers of the dumbbells, 



\ch\ 



\C21 - C43I 



ri + 



0.2 



r2i - rs 



04 



fli + a2 as + 0-4 

as well as the conformation angles 1^21 and 1^4^^, which are given by 



r43 



cos(V^i 



(9) 



(10) 



I ' fij I 

The hydrodynamic interactions between the dumbbells induce a slow rotation of Ch and 
oscillations of D(t) around some mean distance as indicated by figure [2J The frequency 
Q' of these oscillations is the rotation frequency of the dumbbell axes r2i and r^^ with 
respect to Ch- Therefore the frequency of the slow rotation of Ch is equal to the difference 

n-n'. 

b) 





Figure 2. In part (a) the trajectories of the centers of two rotating, equally oriented 
and equally asymmetrically shaped dumbbells are sketched. In addition to that the 
coupling parameters, V21 = V43, are identical. The trajectories of both hydrodynamic 
centers C21 and C43 have the same mean radius, which corresponds to the mean 
length of the vector Ch- Part (b) corresponds to the case of two dumbbells with 
differently asymmetric shapes and different coupling strengths (V21 / V43). The 
resulting trajectories are open and the mean distance D{t) ~ |c/i| increases during 
the motion. 



3.2.1. Two dumbbells having identical shapes and couplings. The dumbbell dynamics 
is considerably simplified for equally asymmetric shapes of the two dumbbells, ai =03, 
and equal couplings, V21 = V43. If the two dumbbells are also equally oriented, as 
sketched in figure |2]^a), the phase angles are equal, 02i = 043 = 0, and the conformation 
angles are identical for all times, il)2i = ip43 = ip- For two oppositely oriented dumbbells 
the phase difference has the constant value A0 = tt whereas the conformation angles 
add up to %Ij2i + '043 = vr. 
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According to the symmetry of the system, the distance D{t) between the 
hydrodynamic centers of two equally oriented dumbbells oscillates around a mean 
distance with the frequency 2Q' as indicated by the sketch in figure EJ^a). The 
corresponding oscillation period is tt/Q', which is shown for a specific parameter set 
in figure [3t^a). In contrast to this, the distance D{t) for two dumbbells with opposite 
shape orientation oscillates with a period of r^' = In being twice as long as in the 
former case. For comparison we show the x component of the driving field /^(t) in figure 
[3l^c), which illustrates the inequality of the frequencies, f2' < fi, and thus the inequality 
of the oscillation periods r^' > tq, = 2t\ j^. 




0.25 0.5 0.75 1 



time / T^, 

Figure 3. Part (a) shows the dimensionless distance D{t)/a for two equally 
asymmetrically shaped and equally oriented dumbbells with ai — — 1.2 and 
V21 = Viz = 0.7. In part (b) we present the periodic contribution to the phase angle 
ip — Tp — Q't. For comparison the periodic variation of the x component of the driving 
field / is shown in part (c). 

In figure El^b) the oscillatory contribution ip to the phase ip = Q't + ip is shown. 
Since the coupling strengths of the dumbbells are constant, the oscillation of about 
a linearly growing part confirms that the rotational hydrodynamic drag changes as a 
function of the relative orientation of the dumbbell axes with respect to c^. In the 
case of equally oriented dumbbells the viscous drag is maximal in the ranges about the 
conformation angles ijj = (2n— l)7r/4 (n = 1, 2, 3, 4) where fluid between the dumbbells is 
squeezed out or sucked into the region between the dumbbells. These are also the phases 
of the rotation during which the repulsion (resp. attraction) between the dumbbells is 
strongest. 

According to the symmetry of the system the attractive forces acting between the 
two dumbbells in one quarter of a period {0 < t < r^r / 4) are compensated during the 
consecutive quarter {tqi/A < t < tqi/2) by repulsive forces of the same magnitude. 
This causes an oscillation of D{t) with the frequency 2Q' around some mean value, 
but without a net attraction or repulsion. In the case of oppositely oriented dumbbells 
two consecutive maxima of the viscous drag have different amplitudes and therefore 



Hydrodynamic attraction and repulsion between asymmetric rotors 



7 



the attractive and repulsive hydrodynamic forces do not cancel each other completely 
during one half of an oscillation period but during a whole rotation with the period r^r. 

During the rotation the viscous drag acting on the inner beads, which are closer 
to the respective other dumbbell, is higher than the drag on the outer beads of the 
dumbbells. These differing drag forces induce an oscillation of the phase angles (pij 
during the rotation. 

In the ranges of increasing (decreasing) viscous drag the rotation frequency Q' of 
the dumbbell axes decreases (increases). When Q' is larger (smaller) than u the phase 
delays (pij of the dumbbells decrease (increase). 

The oscillations of D{t) and (j){t) are shown as a function of time for two different 
values of the coupling constants in figure HI In Part (a) the coupling to the driving field 
is weak, V21 = V43 = 0.7, and in part (b) the coupling is strong, V21 = V43 = 100. The 
amplitudes of the oscillations of the distance are the same in both cases. Is contrast to 
that the phase lag of the dumbbell axes with respect to the orientation of the external 
field is smaller in the case of large values of the coupling. Moreover the diagrams in 
figure m show that the oscillations of are in phase with the oscillations of the distance 
in the case of strong couplings to the external field. Contrary to that (j)(t) is delayed 
with respect to D(t) for weak couplings. 

The reason for this delay is the following: While the distance between the dumbbells 
increases the drag on the dumbbells also grows. So the rotation frequency Q' of the 
dumbbells must decrease. As soon as it is below u the phase angle (j) grows and so 
does the torque on the dumbbells (~ sin0). But for weak couplings the torque builds 
up much slower than for strong couplings. So for small coupling constants even in the 
beginning of the domain where the distance (and thus the drag) decreases (p is still 
growing because the torque is too weak. At some point the torque is large enough to 
overcome the drag so that decreases again afterwards. An analogous argument holds 
in the ranges where D{t) is increasing. The shift between the curves of D{t) and 0(t) 
has to be smaller than a quarter of an oscillation period in our considerations because 
otherwise the field f{t) would be too weak to enforce a synchronous rotation of the 
dumbbells at all. 

For strong couplings there are small dips at the maxima of 0(t). These correspond 
to local minima of the drag, which occur when the dumbbells are aligned with each other 
so that the four beads lie on a single line. In this conformation the distance between the 
dumbbells is maximal. An equivalent phenomenon can be seen in the plot of 0(t) for 
weak coupling. Of course there is also a local minimum of the drag when the dumbbells 
are aligned with each other and the distance is maximal. But due to the weak torque 
acting at that stage there is no dip but only a slight flattening of the curve in this case. 

In the coordinate system rotating with the connection vector the conformation 
after one period tq' maps on top of the original conformation. The reciprocity of 
the motion is illustrated in figure O where the distance D is plotted as a function 
of the conformation angle ip for weak couplings, Vu = V34 = 0.7, in part a) and 
strong couplings, V12 = V34 = 100, in part b). Both curves are invariant under the 
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Figure 4. The distance D{t) (solid line) and the phase lag — (p2i — 043 (dashed line) 
are shown as a function of time for ai — — 1.2. In part (a) the coupling is weak, 
V21 = V43 = 0.7, and in part (b) the coupling constants are large, V21 = V43 = 100. 




45 90 135 180 45 90 135 180 



Figure 5. Trajectories in the D — ip plane for a system of two dumbbells with 
ai — ~ 1.2 and equal coupling constants. In part (a) the couplings are 
V12 = V34 — 0.7 and in part (b) V12 = V34 = 100. Both trajectories are invariant 
under the transformation ip — > —-0 showing the reciprocity of the motion. What is 
more the trajectories fall on top of each other, which illustrates the universality of the 
trajectories. However, according to the two different couplings the rates "0 along the 
trajectories differ. The rates ip are color coded where red denotes a rate below the 
average one and blue a rate above. 



transformation ip — )■ —ip. This underlines the reciprocity of the motion as it can aheady 
be seen in the plot of D{t) which is invariant under time reversal t — ?■ —t. In addition 
to that the curves for weak and strong couplings are identical in the conformation space 
D — Ip as shown in figure [5l In fact the trajectories shown are universal for all values 
of the coupling constants as long as they are large enough so that the dumbbells can 
follow the field f{t) synchronously. 

However, the rate of change of the conformation angle ip along the curves depends 
very much on the coupling parameters (cf. figure H]). This is indicated by the color 
code along the two curves, where red and blue mark the rates ip below and above the 
angular frequency Q'. This shows that irrespective of the strength of the driving field 
the dumbbells pass through the same set of conformations during one cycle, although 
each conformation is passed at a different rate, as indicated by the distribution of the 
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different colors along both lines. The reciprocity of the motion ensures that for a given 
rotational frequency the distance of the dumbbells described by the Stokesian dynamics 
is locked to the individual orientation of the dumbbells with respect to Ch and returns 
to the same value after one period. 




time / time / 



Figure 6. Tlie upper parts show the phase difference A(j) between the two axes of the 
rotating dumbbells as a function of time. The distance D{t) is plotted in the middle 
parts and the lower parts show the x component of the driving field /. In part (a) the 
two dumbbells repel each other. The corresponding parameters are oi = 1.8, 03 = 0.6, 
V21 = 1.2 and V43 — 1.0. In part (b) the dumbbells attract each other. There the 
parameters are ai — 1.2, 03 = 1.1, V21 = 3.2 and V43 — 0.8. 

3.2.2. Two dumbbells of dissimilar shapes or couplings. For two differently shaped 
dumbbells or for different couplings to the external field the phase difference A0 and 
the distance D{t) are plotted as functions of the time in figure |6] for two different sets 
of parameters. For comparison also the x-component fx of the driving field is shown. 
In contrast to the previous section the phase difference A0(t) oscillates in time. The 
time dependence is different for different parameters as indicated in figure Et^a) for the 
parameters ai = 1.8, a^ = 1.0, V21 = 1.2, V43 = 1.0 and in figure [6t^b) for the parameters 
ai = 1.2, as = 1.1, V21 = 1.0, V43 = 2.0. 

In the previous section D{t) was either oscillating with the frequency 2Q' or Q', but 
in both cases D{t) was symmetric with respect to a reflection of time, cf. figure HI Here 
both, D{t) and A0(t), are not symmetric anymore with respect to time reflection. 

As in the case of equally asymmetric dumbbells the viscous drag is a function of 
the angles which the dumbbell axes enclose with the vector Ch. There are also ranges 
of the phase angles, in which liquid is squeezed out between the dumbbells and others 
where fluid is sucked in. As a consequence the length of the vector Ch is also oscillating 
here, i. e. there are phase ranges of dumbbell attraction, which are followed by phases 
of dumbbell repulsion and so on. 

For finite values of A0 7^ 0, which are a consequence of the broken symmetry, 
the conformation angles ip2i and 1^4^^ differ from each other. Combined with the 
hydrodynamic interactions, which depend nonlinearly on the distances between the 




Figure 7. In part (a) the distance D{t) between two dumbbells is plotted as a function 
of the conformation angle '021 between the axis of the first dumbbell and the vector 
Ch- Part (b) shows the derivative of D{t) as a function of 7/121 • The parameters are 
ai = 1.2, 03 = 1.1, V21 = 3.2, and V43 — 0.8. The dynamics of this system is much 
more complex than for equally asymmetric dumbbells with equal couplings (cf. figure 

El). 

le-06 n ' ' ' ' ' n 
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n 

Figure 8. The Fourier modes of the function are shown. There is a non- 

vanishing zeroth order coefficient that leads to an attraction or a repulsion between 
the dumbbells. However, there are eight higher order coefficients that have larger 
values which indicate a quite complex motion. 

beads, the magnitude and the directions of the forces between the beads during one 
rotation cycle are quite complex. Therefore the trajectories in the configurational space 
are not universal anymore. 

For one parameter set the distance D(t) is plotted as a function of the conformation 
angle ^/'2i in figure [Tl^a). It is clear from first sight that the behavior is much more 
complex here than in the case of equally asymmetric dumbbells with equal couplings, 
which was shown in figure [5l Especially the magnitudes of the attraction and repulsion 
between the dumbbells during the different stages of the motion are different here. So 
a priori it is not clear anymore whether the dumbbells can attract or repel each other 
during a whole rotational period. In order to illustrate the complex behavior of D as a 
function of ip2i in more detail, figure [^b) shows the derivative ^^{'ip2i)- In the ranges 
of the conformation angle tp2i where this function is positive the dumbbells repel each 
other and when the function is negative they attract each other. The Fourier modes of 
this function are given in figure [H In this spectrum one can clearly see that there is 
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a non-vanishing coefficient of zeroth order which means that the distance between the 
dumbbells changes during a complete dumbbell rotation. But it is also obvious that 
this net repulsive or attractive effect is superimposed by much stronger oscillations. So 
it is not possible to determine in which phase of the motion the decisive effect takes 
place that causes the overall repulsion or attraction. In fact all coefficients from n = 1 
(~ e*'^^^) up to = 8 (~ e^*'^^^) are larger than the zero order coefficient. 

As shown above the mean distance between two unequally shaped dumbbells may 
increase, as sketched in figure [2]^b), or decreases as a function of time. In figure [9] we have 
plotted D{t) for the same parameter sets as in figure [6l but for a longer time window. 
One can easily see that in the mean the dumbbells repel each other in part a) and 
attract each other in part b). How the hydrodynamically induced repulsive or attractive 
dumbbell motion depends on the parameters, describing the system's asymmetry, is 
illustrated in terms of diagrams for representative sets of parameters in figures [T0lfT2l 




4 8 12 16 4 8 12 16 



time / time / 

Figure 9. The distance D{t) between the centers of the dumbbell is shown as a 
function of the time for the same two parameter sets as in figure [51 In part (a) the 
dumbbells repell each other. The parameters are ai — 1.8, 03 = 0.6, V21 = 1-2 and 
V43 = 1.0. In part (b) the dumbbells attract each other for the parameters ai = 1.2, 
as = 1.1, V21 = 3.2 and V43 = 0.8. 

In figure HU] areas of dumbbell attraction and repulsion are shown as a function of 
the two coupling parameters V21 and V43 for the fixed parameters, ai = 1.2 and = 1.1, 
corresponding to unequal dumbbell shapes. For combinations of the two couplings along 
the solid lines in figure [TOl the bead asymmetries of the dumbbells are compensated by 
the difference between V21 and V43, so that the distance D{t) is constant in the mean. If 
one coupling parameter is much larger than the other one (V21 ^ V43 or V43 ^ V21) the 
dumbbells attract each other. Between those two regimes there is a range of parameters 
for which the dumbbells repel each other. For increasing coupling constants the repulsive 
area widens. 

In figure im the regions of dumbbell attraction and repulsion are shown as functions 
of the two bead radii ai and 03 for fixed coupling constants, V21 = 1.2 and V43 = 1.0. For 
these parameters there is only a narrow domain, in which the dumbbells attract each 
other. In this area the dumbbell asymmetries differ only about 20 — 30%. The phase 
diagram in figure [11] also shows that for identical asymmetries of the two dumbbells, 




1 1.5 2 2.5 3 3.5 



V21 

Figure 10. Ranges of dumbbell attraction and repulsion are shown as a function of the 
coupling constants V21 and V43 for slightly differing dumbbell asymmetries: ai = 1.2, 
03 = 1.1. 

ai = as, one may have either dumbbell repulsion, as in the case ai = = 0.6, or 
attraction, as in the case ai = = 1.8. 




0.6 0.8 1 1.2 1.4 1.6 1.8 



Figure 11. Ranges of dumbbell attraction and repulsion are shown as functions of the 
ratios of the bead radii ai/a2 and a^/a^ for unequal coupling parameters, V21 = 1.2 
and V43 = 1.0. 

In figure fT2] the coupling V43 = 1.0 and the asymmetry as/a^ = 1.1 of one of the 
dumbbells were fixed, whereas the asymmetry ai/a2 as well as the coupling parameter 
V21 is varied. For a symmetric or slightly asymmetric dumbbell, i.e. ai/a2 ~ 1, there 
is always attraction, but for stronger asymmetries there are also regions, where the two 
dumbbells repel each other. Even for equal couplings (V21 = V43) one can find areas of 
dumbbell attraction as well as regions in which the dumbbells repel each other. 

Away from the solid lines of vanishing attraction and repulsion in figure [T0| [TT] or [T2l 
the modulus of the relative velocity between the dumbbells increases with the distance 
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Figure 12. Ranges of dumbbell attraction and repulsion are shown as a function of 
the coupling strength V21 and the bead radius ai. The other parameters are V43 = 1.0 
and as — 1.1. 

from these lines. This is demonstrated in figure [131 where the mean relative velocity 
between the dumbbells is shown as a function of the coupling V43. In this figure the 
same parameter set as in figure [TO] was used and V21 = 1.2 was fixed. In this case the 
transitions from attraction to repulsion and back occur at the values V43 ~ 1.17 and 
V43 ~ 1.59 and there is only a small range of values of V43 for which the dumbbells repel 
each other. 

X 10"* 



-2 

<dD/dt) 
CO a2 

-4 
-6 
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Figure 13. The mean relative velocity between the dumbbells is shown as a function 
of the coupling strength V43 for the same parameter set as in figure [TU] and the fixed 
value V21 — 1.2. 

Is there a simple mechanism by which the attractive or repulsive behavior can be 
explained? As shown before, the broken time reversal symmetry of D{t) is necessary for 
the non-reciprocal motion of the dumbbells to occur. Furthermore for some parameter 
sets the phase shift between the dumbbells as well as its modulation amplitude 
has a minimum in the repulsive regime. However, this depends very much on the 
parameters being used. For other parameter sets A0 takes its maximum value in the 
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range in which the dumbbells repel each other. For a third class of parameter sets 
there is no obvious correlation between the phase shift and the fact that the dumbbells 
behave in an attractive or a repulsive way. All in all the mechanism leading to dumbbell 
attraction or repulsion is a complex interplay between the applied torques resulting from 
the coupling parameters and the drag forces, which depend on the bead radii and thus 
on the asymmetries of the dumbbells. 

4. Conclusions and discussion 

Hydrodynamic interactions between two rotating asymmetric dumbbells in a fluid at 
a low Reynolds number were investigated. We found either a temporally averaged 
attraction or repulsion between the dumbbells if they have different shapes and/or their 
couplings to the driving field differ. Differences in shape and/or coupling break the 
time reversal symmetry, and are hence a pre-condition for attraction or repulsion. No 
generic rule could be identified of whether a specific form of symmetry breaking leads to 
dumbbell attraction or repulsion. We presented phase diagrams separating parameter 
regions of attraction from regions of repulsion. 

We suggest experiments where anisotropic birefringent particles are rotated by 
circularly polarized light, similar to recent experiments [271 [2H], in order to explore 
the hydrodynamic attraction and repulsion between rotating small particles. In such 
experiments the torque on an anisotropic particle may easily be tuned by varying the 
laser power. Different torques may also be applied to the particles by using different 
shapes or materials for the particles. A spatially varying power of the circularly polarized 
light can also cause different torques on neighboring particles. 

Dumbbells or anisotropic particles can also be constructed from super-paramagnetic 
particles which can be rotated using a magnetic field. In this case the magnetic dipole- 
dipole interactions have to be taken into account. Magic angle spinning may suppress the 
magnetic dipole or multipole interactions [25j such that the hydrodynamic interactions 
dominate. 

If the hydrodynamic interactions of an assembly of asymmetric objects are 
investigated, collective dynamics such as chaotic motion can be expected. Our findings 
might find applications in an efficient way to micromix a fluid or to separate particles 
on the small scale. It might also be useful for studies on the collective assembly of 
micro-swimmers. 
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